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Abstract. The effect of redshift on the observation of distant galaxies is briefly 
discussed emphasizing the possible sources of bias in the interpretation of high- z 
data. A general energetic criterion to assess physical self-consistency of evolutionary 
population synthesis models is also proposed, for a more appropriate use of this 
important tool to investigate distinctive properties of primeval galaxies. 
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1. Introduction 
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r- — , , The recent major advances in the observation of the deep Universe 

have substantially improved our chance to reach the redshift of galaxy 
formation. Objects at z > 5 begin in fact to emerge from HST ob- 
servations and other deep surveys with the new-generation telescopes 
(Cohen et al. 2000). 

As far as we move to larger distances, however, optical (and infrared) 
observations probe galaxy spectral energy distribution (SED) at much 
shorter restframe wavelength, in the ultraviolet domain. If neglected, 
this change in our spectral vantage point could lead to an important 
bias in the interpretation of high-redshift data. In this note, we will 
briefly analyze some aspects of this problem dealing both with the 
observational side and the appropriate use of theoretical tools such as 



population synthesis models to match primeval galaxy evolution. 



2. Redshift constraints and galaxy spectral properties 

One important consequence of redshift (besides the first and most ob- 
vious effect of moving spectral features to longer wavelength) is that 
the wavelength interval AX h s sampled by any fixed set of photomet- 
ric bands will span a narrower portion of the restframe SED when 
observing galaxies at increasing z, such as AX res t = AA & S /(1 + z). 
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Figure 1. Restframe portion of galaxy SED explored at different redshift by the 
Johnson four-color photometry through the B, V, R, and I bands. The wide wave- 
length range, between 4000 and 9000 A, attained for local galaxies, "shrinks" to 
a 1000-2200 A interval when observing objects at z — 3. The synthetic SED of a 
model elliptical is displayed, for reference. 



This effect, by itself, can be of paramount importance as far as we 
try to gain information on galaxy stellar content by means of multicolor 
photometry. Figure 1 gives an illuminating example in this sense. 

For a z = galaxy, for instance, Johnson four-color photometry in 
the BVRI bands would span a wide wavelength range, from 4000 to 
9000 A, and collect a substantial contribution from the galaxy stellar 
population. On the contrary, when observing at z = 3 with the same 
photometric system, we would be covering a scarce ±600 A spectral 
window centered about 1600 A in the galaxy restframe. Quite deceiv- 
ingly, our "multicolor" photometry at z = would translate into a 
nearly "monochromatic" estimate of the galaxy flux at z = 3. 

2.1. Sampling stellar populations of high-z galaxies 

According to the overall color-magnitude distribution, each star in a 
galaxy aggregate contributes with a different weight to the integrated 
luminosity with varying wavelength. In Buzzoni (1993) we approached 
this problem in a statistical way introducing the concept of "effec- 
tive" contributors in a stellar population. This quantity directly relates 
to galaxy integrated properties, such as total luminosity or surface 
brightness fluctuations (as further explored, for instance, by Tonry and 
Schneider 1988). 

We showed, in particular, that by sampling at a given wavelength 
a total luminosity Ltot provided by N stars, the expected Poissonian 
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Figure 2. Relative number of "effective" stellar contributors {N e ff) and SED (Ltot) 
for a 5 Gyr SSP model with Salpeter IMF and [Fe/H] = +0.3 dex according to 
Buzzoni (1993). Both quantities are in logarithm scale and arbitrarily normalized 
to their values at 1000 A. 



fluctuation results: 
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where li is the luminosity of each composing star. 

If all stars have the same luminosity, then simply a(L)/Lt t 
l/yN. More generally, however, a quantity 



N 



N P . 



(E JV U) 
E 



n 1? 



< N 



(2) 



can be defined, as a function of A. Its value gives a statistical measure 
of the effective number of stars that provide L tot at a given wavelength. 
Following Buzzoni (1993), a "mean" luminosity (l e ff) for the composing 
stars can also be computed as 
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so that, at every wavelength, 

Ltot = N eff x l t 



(3) 



(4) 



In Fig. 2 we show the expected value of N e ff for the illustrative 
case of a 5 Gyr simple stellar population (SSP) with Salpeter IMF and 
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[Fe/H] = +0.3 along with the integrated SED of the model. This should 
closely resemble the real case of primeval elliptical galaxies (Buzzoni 
1995). 

As a stricking feature of the model, note that below 2000 A the 
effective number of "luminous" contributors dramatically drops by over 
six (!) orders of magnitude compared for example with the B band 
value at 4000 A. On the contrary, luminosity over the same wavelength 
range only dims by a factor of Lb /L2000 ~ 10 3 . In terms of surface 
brightness fluctuations this means, for instance, that the ultraviolet 
image of a galaxy consisting of 10 11 stars should fairly compare with 
the B observations of a (10 6 star-poorer) Galactic globular cluster. 

Such a severe undersampling of the galaxy stellar population in 
the ultraviolet range could have an important impact when studying 
morphology of z > 1 ellipticals at optical wavelength (but the same 
argument holds also for late-type systems). In particular, high-redshift 
objects would more likely appear as coarser systems with increasing z 
(e.g. Puerari 2001, this conference). 



2.2. Redshift and age bias 

The selective sampling of galaxy stellar population with changing wave- 
length has even more pervasive consequences, as far as we track evolu- 
tion of star forming systems at different redshift. 

More generally, we could describe a composite stellar population in 
terms of a suitable convolution of SSPs according to the star formation 
rate (SFR) at the different epochs. The total luminosity at age "£" 
results therefore: 

Ltot{t) = [ L SS p(t) SFR(t - r) dr (5) 

Jo 

A mean luminosity-weighted age of the composing SSPs can also be 
derived such as 

f*TL SSP (T)SFR(t-T)dT 

u = « (6) 

This could be regarded as the "representative" age of the stars that 
contribute to galaxy luminosity. We know (Buzzoni 1995; Tinsley and 
Gunn 1976) that SSP luminosity evolution can be pretty well described 
by a simple power law such as L$sp °c t~ a with the power index a that 
is a function of wavelength. As an instructive example, if we consider for 
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the SFR a power-law time decay such as SFR oc t"^, 1 eq. (6) becomes 

jt T l-a {t _ T) - VdT 



' / *r-«(t-r)-"dr 
The integral has an analytical solution such as 



(7) 



- 1 — a 

U = - t. 8) 

2 — a — 77 

In bolometric, a = 0.76 (Buzzoni 1989) and, for a constant SFR (i.e. 
Tj = 0), t* = 0.2 t. At 15 Gyr, bright stars are therefore in average 3 
Gyr old, but this value changes according to the a(\) function. For a 
SSP with Salpeter IMF and Zq, the luminosity power index smoothly 
increases from a = 0.71 in the K band to 0.94 in the B band, and 
exceeds unity below 3500 A. A value of a > 1 implies that i* — ► 0, that 
is only the youngest and more massive stars are representative con- 
tributors to galaxy ultraviolet luminosity, as discussed in the previous 
section. 

The mean age of stars contributing to galaxy luminosity at different 
wavelength for a synthesis model of a 15 Gyr Magellanic irregular (see 
Buzzoni 2001 for details) is reported in Fig. 3. Note from the figure 
that the value of i* smoothly decreases with decreasing wavelength. As 
a result, high-redshift objects forcedly appear to be younger than local 
homologues at z = in spite of any intrinsic evolution. 



3. Energetic self-consistency of population synthesis models 

Evolutionary population synthesis models are extensively used as a 
reference tool to reproduce galaxy colors and investigate galaxy spec- 
trophotometric properties at the different cosmic epochs (Buzzoni 1989, 
1995, 2001; Bruzual and Chariot 1993; Worthey 1994; Bressan et al. 
1994; Fioc and Rocca-Volmerange 1997). 

Besides any difference in the computational details and input physics, 
a crucial constraint that needs to be consistently fulfilled in the syn- 
thesis procedure concerns a suitable match of the MS and Post-MS 
luminosity contribution in the SSP models. The so-called "Fuel con- 
sumption theorem" of Renzini and Buzzoni (1986) provides a simple 

1 The major advantage of this parameterization is that the full range of galaxy 
morphologies can simply be accounted for by an age-independent distinctive 
birthrate: b = SFR/ < SFR >= (1 — 77). This easily derives recalling that, at 
age t, SFR = C t~ v and its time average is < SFR >= C t^ 1 f t^ 71 dr. Providing 
that r\ < 1, we always have < SFR >= SFR/(1 — 77), from which the value of b 
directly follows, by definition. 
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Figure 3. Theoretical SED [upper panel), and "mean" age (relative to galaxy age) 
of stars at different photometric bands (and in bolometric as well), according to 
eq. (7) (lower panel), for a 15 Gyr Im model galaxy from Buzzoni (2001). 



and very powerful tool in this sense. According to the theorem, the 
Post-MS bolometric luminosity in a SSP of total luminosity (Lb {) is 



L 



PMS 



B x L hn i x Fuel. 



(9) 



In the equation, the nuclear fuel refers to the amount consumed by 
one star leaving the MS turn-off (TO) point along the whole Post-MS 
evolution. This is a natural output of the theoretical stellar tracks and is 
therefore known with high accuracy. The scaling factor B is the "specific 
evolutionary flux" of the SSP, and directly deals with the stellar clock, 
that is the relationship between MS lifetime and TO stellar mass (see 
Renzini and Buzzoni 1986). In particular, we have that 



B oc MZ% x M T o 



(10) 



assuming that a power-law IMF holds with dN* oc M~ s dM* (s = 2.35 
for the Salpeter case). Again, the r.h. term of eq. (10) is a direct output 
of stellar evolution theory, and B should therefore be regarded as a 
"universal" function, that is nearly insensitive to the specific input 
physics of the SSP models. 

A comparison of the B scaling factor, as derived from the synthesis 
codes of Bruzual and Chariot (1993), Buzzoni (1989), and the Padova 
group (Bressan et al. 1994) for a SSP with Zq and Salpeter IMF is 
displayed in Fig. 4. While both the Buzzoni and Padova codes consis- 
tently match the fuel consumption theorem, the Bruzual and Chariot 
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Figure 4- The "specific evolutionary flux", B, for SSP models with Zq and Salpeter 
IMF according to Buzzoni (1989; "•" markers), Bressan et al. (1994; "□") and 
Bruzual and Chariot (1993; "o"). The value of B is given in unit of yr -1 Lg 1 . It 
directly scales Post-MS contribution to SSP luminosity through eq. (9). 



models show a much higher value for the specific evolutionary flux. 
The effect is a macroscopic one, and reaches about a 100% discrepancy 
at intermediate age. Such higher scale factor leads, in eq. (9), to a 
corresponding overestimate of Post-MS luminosity at every age. The 
effect is quite shifty as it will not explicitly appear when considering 
luminosity evolution (cf. Fig. 5, left panel), but would strongly affect 
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Figure 5. Left panel: bolometric evolution for a SSP with Zq and Salpeter IMF 
according to different synthesis codes, as labelled. Luminosity is normalized at 12 
Gyr. Right panel: apparent B — V evolution of elliptical galaxies according to Buzzoni 
(1995; "Buz") and Bruzual and Chariot (1993; "B&C"). A cosmological model with 
(H , q , Zf) — (50 km s _1 Mpc -1 , 0, 30) is assumed. The absolute age of galaxies at 
the different redshift is also labelled on the plot. 
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the theoretical prediction of galaxy M/L ratio and apparent color evo- 
lution. In the latter case, enhanced Post-MS contribution would lead to 
a larger number of red giant stars and therefore a systematically redder 
color for the stellar aggregate as a whole. This is shown in Fig. 5 (right 
panel) when comparing the Bruzual and Chariot (1993) results with the 
expected evolution of elliptical galaxies according to Buzzoni (1995). 
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